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Abstract. A basic assumption of tiling theory is that adjacent tiles can meet in only 
a finite number of ways, up to rigid motions. However, there are many interesting tiling 
spaces that do not have this property. They have "fault lines" , along which tiles can 
slide past one another. We investigate the topology of a certain class of tiling spaces of 
this type. We show that they can be written as inverse limits of CW complexes, and 
their Cech cohomology is related to properties of the fault lines. 



I. Introduction 

In discussing tilings, a standard assumption is that tiles can meet in only a finite 
number of ways, up to rigid motion. Equivalently, for any radius R, there are only a 
finite number of patches of radius R in the tiling, up to rigid motion. This condition is 
called finite local complexity, or FLC. Q 

In this paper we consider tilings that do not meet the FLC condition. We show that 
spaces of such tilings can be given a natural topology in which they are compact. Many of 
the techniques used for FLC tilings, such as inverse limit constructions and cohomology 
calculations, can be modified to handle non-FLC tilings. In particular, we work out a 
number of simple examples, and prove theorems about a broader class of examples. 

1.1. Past work. The first substantial work on non-FLC tilings was done by Kenyon 

(a a a' 
a a a 
a a a, 

but in which one of the columns was shifted by an irrational distance relative to the other 
two. (The tile itself was not a polygon. Rather, it had straight edges on the left and 
right, while the top and bottom had fractal pieces that looked like the devil's staircase). 
Upon further substitution, the shift between the columns would grow to an arbitrarily 
large multiple of the original irrational distance. In the limit, an infinite line of tile edges 
would appear, along which tiles could face one another in an infinite variety of ways. 
Such a line is an example of a "fault line" , which we now define. 

Given a finite set of prototiles (usually assumed to be polygonal, or perhaps closed 
topological disks), a tiling is a covering of M d by rigid motions of copies of these prototiles 
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^Frequently an even stronger condition is applied, namely that tiles can only meet in a finite number 
of ways up to translation, a condition that excludes tilings like the pinwheel [Rad , in which tiles appear 
in an infinite number of orientations. 

1 



2 



Frank and Sadun 



that are only allowed to intersect on their boundaries. A tiling space X is any translation- 
invariant set of such tilings that is closed under the third topology of section 11.21 Let 
T G X be a tiling containing an infinite line (or ray) i of tile edges, and let v be a unit 
vector parallel to i. We say that t is a fault line if for every e > there is a tiling T'el 
such that: 

(1) On one half- (or quarter-) plane with boundary £, T' = T, and 

(2) there is a t with < \t\ < e such that on the other half- (or quarter-) plane with 
boundary i we have T = T' + tv. 

If for every sufficiently small e > it is possible to choose t = e, we call i a regular fault 
line. 

Fault lines play a central role in any discussion of non-FLC tilings. This is because 
Kenyon also proves |Ke2j that if a tiling that is made from a finite prototile set has an 
infinite number of inequivalent two-tile patches, then those patches occur along straight 
edges, or they occur along an entire circle of tile boundaries. The former case leads 
to fault lines, while the latter can only occur for special prototile sets (and never, for 
instance, in primitive substitution tilings [FRj). 

In 1998, Sadun [SaJ proposed some generalizations of the pinwheel tiling. One such 
example (Til (1/2)) uses polygonal tiles (two similar right triangles) that appear in in- 
finitely many orientations, and has a fault line. Interestingly, this example meets the 
conditions of Goodman-Strauss' matching rules theorem [GS] . To wit: there exists a 
finite collection of tiles and a finite set of local matching rules such that these tiles tile 
the plane, but only in a manner that is locally equivalent to the generalized pinwheel. 
That is, a finite set of local rules forces a global hierarchical structure that in turn forces 
infinite local complexity. 

Danzer |Daj extended the theory of FLC tilings and provided additional examples. 
Finally, Frank and Robinson [22 IFRj have considered a large family of "direct product 
variation" (DPV) tilings. These are obtained from products of 1-dimensional substitu- 
tions by rearranging the positions of the tiles within an order-1 supertile. The examples 
of this paper are all DPV tilings in which the rearrangements are purely horizontal, 
thereby preserving the decomposition of the tiling into rows. 

In all cases, the essential feature that prevents FLC is the presence of a fault line. 
Along the fault line, the evolution is described by two 1-dimensional substitutions. One 
describes what happens along one side of the fault line, while the other describes what 
happens along the other side. The two substitutions cannot be Pisot with the same 
stretching factor, because in that case the FLC property is preserved |FRj and there is 
no possibility of a fault line. However, if the substitutions are not Pisot, then differences 
in distributions of lengths on opposite sides of the fault line will generally grow with 
successive substitution, and the FLC property can be lost. This will be discussed further 
in section [2J 

1.2. Three Topologies. There are three metrics, and hence three topologies, that are 
frequently applied to tilings and tiling spaces. 

In the first metric, two tilings are considered e-close if they agree, up to a translation 
of size e or less, on a ball or radius 1/e around the origin. In this topology, the closure 
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Figure 2. A patch of a DPV tiling with vertical and horizontal fault lines. 



4 



Frank and Sadun 



of the (translational) orbit of a tiling is compact if and only if two conditions are met: 
(a) there are only a finite number of tile types, up to translation, and (b) tiles can only 
meet in a finite number of ways, up to translation. This is the most frequently used 
topology for tilings that meet these conditions. Obviously, tiling spaces without FLC 
are not compact in this topology, but neither are pinwheel-type spaces, which may have 
FLC even though the tiles may appear in infinitely many orientations. 

In the second metric, one applies a metric to the group of rigid motions of the plane 
(e.g., defining an e motion to be a translation of size < e followed by a rotation by < e 
about the origin). One then considers two tilings to be e-close if they agree, up to an 
e motion, on a ball of size 1/e about the origin. For tilings in which tiles appear in 
only finitely many orientations (e.g., the Penrose tiling), this is equivalent to the first 
topology. However, it also handles pinwheel-like spaces appropriately. In this topology, 
the orbit closure of a tiling is compact if and only if the tiling has FLC. 

Finally, one can use a metric in which two tilings are e-close if they contain the same 
tiles out to distance 1/e, and if each tile in the first tiling is within an e-motion of the 
corresponding tile in the second tiling. Note that in this topology, an e shear along a 
fault line yields a tiling that is e-close to the original, while in the first two topologies it 
does not. 

For tilings with FLC, the third topology is the same as the second, insofar as it is 
impossible to apply a small rigid motion to one tile without applying the same motion 
to all of its neighbors. However, in the third topology the orbit closure of any tiling 
with finitely many tile types is always compact. (To see sequential compactness, start 
with an arbitrary sequence of tilings, pick a subsequence in which the type, location and 
orientation of a tile near the origin converges. Then pick a subsequence in which the 
type, location and orientation of a second tile converges. Keep working outwards from 
the origin, and then apply a Cantor diagonalization argument to find a subsequence that 
converges everywhere.) 

In the third topology, FLC is not a topologically invariant property. Radin and Sadun 
[RSJ constructed a pair of spaces, one FLC and one not, that are topologically conjugate. 

1.3. Outline of paper. In section[2J we study the evolution of fault lines in 2-dimensional 
substitution tilings. This is essentially 1-dimensional dynamics, and we relate properties 
of the fault line to the form of the induced 1-dimensional substitutions. 

In section [3l we consider a 2-dimensional substitution tiling with horizontal fault lines. 
The rows of this tiling are (almost) all the same, up to a horizontal shift, which is 
controlled by a vertical 1-dimensional substitution. We show that the resulting tiling 
space can be constructed as the inverse limit of compact CW complexes. We explicitly 
compute the cohomology of this tiling space, and discuss the meaning of each term. 

In section |U we consider a more complicated substitution, as a step towards the direct 
product variations considered in section These direct product variations look like 
the product of a vertical and a horizontal 1-d substitution tiling, except that the rows 
are sheared by an amount governed by the vertical substitution, and exhibit horizontal 
fault lines. The cohomology of the resulting tiling space is computable in terms of the 
cohomologies of the vertical and horizontal 1-d spaces, and the combinatorics of the 
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vertical substitution. Specifically, let /i be H l of the horizontal 1-d substitution tiling 
space, and let M be the n x n substitution matrix of the vertical substitution, as applied 
to collared tiles. Then 

• H° of the 2-d tiling space is Z, of course. 

• H 1 of the 2-d tiling space is isomorphic to H 1 of the vertical 1-d substitution 
space. 

• H 2 of the 2-d tiling space is isomorphic to the tensor product of [i with the direct 
limit of Z n under the map M T . This in turn is related to H 1 of the vertical 
substitution space and the number of possible fault lines. 

• H 3 contains one copy of fi £g> fi for each possible infinite fault line. 

• H k is trivial for k > 3. 

Since if 3 is nontrivial, the 2-d tiling space is not homeomorphic to any 2-d tiling space 
with FLC. 

Finally, in section O we consider open problems in the theory of tilings without FLC, 
and discuss our partial understanding of these problems. 



2. Analyzing a fault line — 1 dimensional dynamics 
Consider the 1-dimensional substitutions <J\(a) = ba, <J\{b) = aaa, 02(a) = ab, 
02(0) = aaa. Both substitutions have substitution matrix ( ^ ^ j , with Perron- Frobenius 

eigenvalue A = (1 + vT3)/2 « 2.3028, the larger root of the equation A 2 — A — 3 = 0. For 
a self-similar tiling, the a tile can be given length A while the b tile has length 3. Note 
that for any word W, a^W) is a cyclic permutation of ai(W), obtained by removing an 
a from the end and sticking it on the beginning. If W is a bi-infinite word, then a\{W) 
and 02 {W) are the same, up to translation by the length of a. This implies that the tiling 
spaces defined by 01 and 02 are exactly the same. 

Suppose we have a tiling with rectangular tiles a and b of widths A and 3, respectively, 
and suppose that cr 1 acts as a substitution on lower edges and 2 acts as a substitution 
on upper edges. Let us construct a horizontal fault line, where the evolution above the 
line is governed by <Ji and the evolution below the line is governed by 02. If at some 
stage there is a pair of exactly aligned a tiles, one above the line and one below, then on 
successive substitutions we will see 



/..s fba\ ( aaaba\ fbababaaaaba\ ( aaabaaaabaaaababababaaaaba 
\ab J ' \abaaa J ' \abaaaababab J ' \abaaaababababaaaabaaaabaaa 

Note that in the first and third substitution, the a tiles are found more on the right of 
the top row and on the left of the bottom row, while in the second and fourth substitu- 
tions they are found more on the left of the top row and the right of the bottom row. 
The reason is that the difference between the number of a tiles up to a certain point 
grows like the second eigenvalue of the substitution matrix, namely 1 — A ~ —1.3. As 
we continue to iterate, this discrepancy grows without bound. (Strictly speaking, the 
discrepancy gets multiplied by 1 — A each time and then adjusted by 0(1) edge effects. 
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Once the discrepancy grows beyond a certain point, the edge effects are dominated by 
the multiplicative factor of 1 — A and we have exponential growth in the discrepancy as 
a function of the number of substitutions. ) 

Pick a point along the fault line. If there are m more a tiles in the top row than the 
bottom up to that point, then the left edges of the tiles on the top row will be offset by 
Am (mod 3) relative to the left edges of the tiles on the bottom row. By continuity, the 
discrepancy takes on all integer values between and m as we move from the left edge 
of the pattern to the point in question. Since m is unbounded and A is irrational, this 
means that the possible offsets of tiles in the top and bottom rows takes on a dense set 
of values in the limit of infinite substitution. In fact the left endpoints of upper a tiles 
are dense in the lower b tiles, because the only way for the discrepancy to grow from m 
to m + 1 is for an additional a tile to appear along the top with a b tile below it. Thus 
every possible adjacency between an upper a an a lower b can occur in the orbit closure; 
by primitivity of the substitutions this implies that any adjacency is possible between 
any upper and lower tiles, this, then Thus we have not just a fault line, but a regular 
fault line. 

Note how the form of the fault line depends on the second eigenvalue of the substitution 
matrix. If the substitution were Pisot, then the discrepancy in the number of any species 
of tile would be bounded, and the offsets between tiles would take on only a finite number 
of values. As a result, the FLC condition would be preserved. For a detailed proof that 
Pisot substitutions do not lead to fault lines, see [FR] . 

Finally, note that there are only two possibilities involving substitutions on two letters. 
If the discrepancy in the number of a tiles grows without bound, then we have a regular 
fault line. If the discrepancy is bounded, then we preserve FLC. It is not known whether 
irregular fault lines, in which the set of possible offsets is infinite but not dense, are 
possible in substitution tilings. In any case, they would require more than two letters. 



3. A SIMPLE 2-DIMENSIONAL EXAMPLE 

Consider a 2-dimensional tiling with two rectangular tiles. Both the A and B tiles 
have height 1, but the A tile has width A = (1 + vT3)/2 and the B tile has width 3. 



We consider the self-affine substitution Ti(A) = 




any n e N, an n-supertile is a collection of tiles of the form £ n (v4) or E n (S). A tiling 
is allowed by the substitution if each of its finite patches of tiles can be found in some 
n-supertile. The smallest closed set (under the third topology) containing all allowed 
tilings is called the substitution tiling space X^,. We will see that measure-theoretically, 
almost all of the tilings in are allowed by the substitution, but that the ones that are 
not are the ones that make the topology different than in the FLC case. 

Note that whenever two supertiles meet along a horizontal boundary, applying £ 
changes the bottom of the top supertile by o\ and the top of the bottom supertile by 
cr 2 - By the results of section [2, the substitution tiling space defined by £ exhibits 
horizontal regular fault lines. 
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Not every row is subject to arbitrary shears. The rows themselves are labeled by points 
in the dyadic integers, describing their hierarchy in the vertical substitution. The label 
in the nth spot is if the row is in the lower (n — l)-supertile of its n-supertile and 1 if it 
is in the upper (n — l)-supertile. If the labels of two adjacent rows differ only in the first 
digit, then the dyadic label of the upper row begins with a and the label of the lower 
row begins with a 1. One can see that the sequence of tiles in the two rows are identical, 
with the upper row offset horizontally by A. If they differ only in the first two digits, 
then the dyadic label of the upper row begins with 01 and the label of the lower begins 
with 10, and the upper row is the same as the lower row, but offset by A 2 — A. If they 
differ only in the first three digits, the upper and lower labels begin with 001 and 110 
resp., and then the rows are offset by A 3 — A 2 — A. If they differ in the first n digits (and 
agree thereafter), then the upper and lower labels begin with n_1 l and 1™ _1 resp., and 
they are offset by A™ — A n_1 — • • • — A. (In general, one sees the new offset as A times 
the previous offset, minus A.) However, in some tilings there exists a row with dyadic 

label 1111 .. . and an adjacent row above it with label 0000 These rows do not have 

to have the same sequence of tiles, and their offset is arbitrary. 

Put another way, all tilings in the tiling space contain horizontal lines separating iden- 
tical rows of tiles, offset by arbitrarily large amounts. However, in a small set of tilings 
(corresponding to a single orbit in the dyadic solenoid that describes the vertical hierar- 
chical structure) there exists a fault line in which the tiling above the fault is unrelated 
to the tiling below the fault. These special tilings have measure zero with respect to all 
translation-invariant measures, and hence have no effect on measure-theoretic properties 
of the tiling space. 

However, they have a tremendous effect on the topology of the tiling space. Indeed, 
this part of the tiling space is 3-dimensional! To describe such a tiling we must specify 
the height of the infinite fault line (a point in M 1 ), the row immediately above the fault 
line (a point in a 1- dimensional tiling space) and the row immediately below the fault 
line (a point in a 1-dimensional tiling space). 

3.1. The approximant L. We will show how is the inverse limit of an approximant 
L under a bonding map induced by X (which we will again call E). The CW complex L 
is actually 4-dimensional , but S : L — > L is not onto, and S(L) is only a 3-dimensional 
subset of L. We will see that the inverse limit of L under E exhibits the right combination 
of 2- and 3-dimensional elements. 

Let o = o"i, and let K be an Anderson-Putnam complex of the 1-dimensional substi- 
tution a, obtained by using collared tiles with a sufficiently large radius D. In order to 
ensure that a 2 is a shift of o"i, we pick D > (A + 3)/(A — 1) so that \D > D + A + 3. The 
1-dimensional tiling space X a is the inverse limit of K under a bonding map induced by 
a. We then let 

(2) L = KxKxKx [0,1]/ ~, (x, y, z, 0) ~ (w, x, y, 1). 

This is understood as follows. Of the coordinates (x, y, z, t) of a point in L, t describes 
the height of the origin in the row containing the origin, and runs from to 1. The 
variable y G K describes a horizontal neighborhood of size D around the origin. In 
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Figure 3. The data encoded in the approximant L. The variables x, y 
and z are points in K that describe the rows nearest the origin, while t is 
the height of the origin in its row. 



other words, it describes the row containing the origin. The variables x and z similarly 
describe the rows immediately below and above, respectively. If the origin sits exactly 
on the boundary of two rows, we may describe the situation either as (x, y,z, 0), with 
the origin sitting on the bottom of the "y" row, or as (w,x,y, 1), with the origin sitting 
on the top of the 11 x" row. Under the identification ~, information about w or z is lost, 
and we know only about the two rows touching the origin ("x" and "y"). 
The bonding map induced from the substitution E is 



(3) H{x,y,z,t) 



{a 2 (y),a 1 {y),a 2 (z),2t-l) ] if t > 1/2, 
(a 1 (x),a 2 (y),a 1 (y),2t); if t < 1/2. 



However, <r 2 (y) is just a translate of <7i(y) = cr(y), and likewise for <r 2 (z), so we can 
rewrite this as 



(4) X(x,y,z,t) 



(a(y) + A, <i(y), a(z) + X,2t — 1); if t > 1/2, 
( ( r(x), ( 7(y) + A, ( 7(y),2t); if t < 1/2. 



Note that, depending on the value of t, the information from either x or z is lost, so that 
the image of S is 3-dimensional. 

Every translate is homotopic to the identity map, so the map E is homotopic to 



(5) T,'(x, y, z, t) = 



(<j(y),<j(y),<j(z),2t-l); if t > 1/2, 
(a(x),a(y),a(y),2t); if t < 1/2. 
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3.2. Computing H*(X^). To compute the cohomology of Xs, we simply compute the 
direct limit of H*(L) under £* = (£')*. Unlike E, E' factors as a product of a vertical 
and horizontal map: E' = £1 o £ 2 = E 2 o E 1; where 

(6) Ei (a;, y, z, t) = (a(x), a(y), a(z),t), 

w ^ ((y,y,z,2t-l); if t > 1/2, 

(7) S 2 {x,y,z,t) = < -f+^i/o 

[(x,y,y,2£); if £ < 1/2. 

Since Si and E 2 commute, the direct limit of H*(L) under E* = EJ o £;?; can be 
computed in two steps. First we take the direct limit under £*, and then we take the 
direct limit under E 2 . Let ft = if 1 (if), and let /i = H l (X a ) be the direct limit of ft under 
a*. Our strategy is 

(1) Using Mayer- Vietoris, compute H*(L) in terms of ft. 

(2) Take the direct limit of H*(L) under E^ Since Ei is essentially just the horizontal 
substitution a, this merely replaces each occurence of ft with [i. Note that we 
never have to explicitly construct K or compute ft. 

(3) Finally, take the direct limit under E 2 . 

Step 1. We take V to be a neighborhood of £ = (say, the set t < 0.2 U t > 0.8) 
and U to be the region where t is not close to zero (say, 0.1 < t < 0.9). U retracts to 
K x K x K x {0.5}. Let ft x be the pullback of ft from the x factor, and likewise for ft y 
and ft z . We then have 

H°(U) = Z, 

H\U) = ft X ® fty® ftz, 

(8) H 2 (U) = {ft X ®fty)@{ft X ®ft Z )®{fty®ft Z ), 

H 3 (U) = ft X ®fty®ft Z . 

Likewise, V retracts to K x K x {0} ~ K x K x {1}. Let ft xy and ft yz denote the pullback 
of ft from the first and second factors, respectively. That be viewed either as 

ft x from t = or ft y from t — 1. We then have 

(9) H°(V) = Z, H\V)=ft xy ®ft yZ) H\V)=ft xy ®ft yZ) H 3 (V) = 0. 

The intersection U D V retracts to two copies of K x K x K, say one at £ = 0.15 and at 
£ = 0.85, and we have 

H°(unv) = z 2 , 

H\Ur\V) = (ftx® fty® ftz)® (ftx® fty® ftz), 

(10) H 2 (UDV) = ((ftx ® fty) ® (ftx ® ftz) ® (fty ® ftz)) 

®((ftx ® fty) ® (ftx ® ftz) ® (fty ® ftz)) 
H 3 (U<1V) = (ftx® fty ® ftz) ® (ftx® fty ® ftz)- 

The Mayer- Vietoris sequence is 

(11) ► H k (L) -A H k (U) ® H k (V) /J fe (f/ nV)^ i/ fc+1 (L) 
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where p is restriction and v is signed restriction. Using a basis for p to make bases for 
H k (U), H k (V) and H k (U PI V), and writing the lower copy of the basis of H k (U PI V) 
first, the matrices of v are 



(12) 



(13) 







1 






V - 


■( 
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(14) 



(15) 
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on H 2 , 



-1 / 



on H 3 , 



where the l's are actually identity matrices that depend on the dimension of p. Note 
that these are all injective, except in dimension 0. As a result, the maps p must all be 
zero (except in dimension 0), and for k > we have that H k (L) is the cokernel of the 
previous v. To summarize, 



(16) 



H°(L) 
H\L) 
H\L) 
H\L) 
H\L) 



Z. 

z, 

(fi <g> p) © (ft <g> p), 
p® p® p. 



By computing d* of the generators of H k (U D V), we see that the generators of H*(L) 
are 1 in dimension 0, dt in dimension 1, dx U dt = dy U dt = dz U dt in dimension 2, 
dx U dy U dt = dy U dz U dt and di U (iz U (it in dimension 3 and ofo; U dy U U dt 
in dimension 4. Here we have used dx, dy, dz and dt as shorthand for 1-dimensional 
cohomology generators in the p x , p y , p z and circle directions. 

Step 2. Taking the direct limit under merely converts p to p. Using 1-dimensional 

methods [BD2J p is easily shown to be the direct limit of Z 2 under the matrix 



and is isomorphic to Z[l/A]. 
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Step 3. From equation (J2j) we compute the effect of £2 011 our cohomology generators. 

= 1 

S*(dt) = 2dt 

(17) Z* 2 (dyUdt) = 2dyUdt 

£2 ((&£ U dy U c?t) = dxUdyUdt 

£2 (da; UdzU c?t) = <ix U % U (it + c/y U cfe U (it = 2<ir U cfa/ U (it, 

E^dz U dy U U dt) = 0, 

which gives us 

#°(X S ) = Z, 

^(Xs) = Z[l/2], 

(18) # 2 (X S ) = p®Z[l/2] = //[l/2], 

# 3 (X S ) = ^/i, 

# fc (X s ) = for k > 3. 

Note that the identifications at the fault line prevent there being any contribution of fx 
to H 1 (Xy 1 ). We do get contributions from ft to H 2 and H 3 . In particular, the H 3 term is 
easy to understand. One factor of fi comes from the tiling above the fault line, one factor 
of fx comes from the tiling below the fault line, and the dt term comes from the location 
of the fault line. In this example, H 2 (Xy) equals the tensor product of fx with H l of the 
dyadic solenoid; however, this is not a general pattern. In the next section we shall see 
an example in which H 2 of the 2-dimensional tiling space is not the tensor product of fx 
with H l of the vertical substitution space, and in section \5\ we will compute a general 
formula for H 2 (X%). 

Finally, note that the form of the answer had nothing to do with the details of the 
substitution a, except that its expansion constant is not a Pisot or Salem number. Let Wi 
and W2 be any two words in the letters a and b, and consider the substitutions 61(a) = w\a, 
62(b) = W2d, 62(a) = awi, 62(b) = aw2- The substitutions 6\ and 62 generate the same 
1-dimensional tiling space. If the pair (61, 62) generates a fault line, as in section [21 then 

consider the 2-dimensional substitution Q(A) = ( ^ ^} J, 0(B) = ( ^ 2 V where 



Wi A J' ^ v 1 \W 2 A ^ 
W\ and W2 are the same as W\ and W2, only written in capital letters. 6 gives rise to a 
2-dimensional tiling space with horizontal fault lines, and the calculation of this section 
can be repeated, line by line, to show that the cohomology of Xq is identical to that of 
Xs, only with fx replaced by ^(Xq). 



4. A MORE INVOLVED EXAMPLE 

Our next example is a direct product variation, in the sense of Frank [Ff] . The vertical 
factor is the period-doubling substitution — > 01, 1 — > 00, while the horizontal factor 
is our usual 1-dimensional substitution a. We let A and B denote a ® 1 and b ® 1, and 
abuse notation by referring to a <g> and b ® as a and b, respectively. The result is a 
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Figure 4. The Anderson-Putnam complex for the period-doubling substitution 

tiling with four rectangular tiles. The tiles A and a have height 1 and width A, while B 

and b have height 1 and width 3. Our substitution is 

(19) 

E(a)=(i B ), E(A)=(l b ), m=( A A A ), nB)=( a a a )- 

v ' yb a J v / yb a J w ya a a J v / ya a a J 

Note that if we ignore the difference between capital and lower case letters, we revert 
to the example of section [3j The period-doubling substitution space is an almost 1- 
1 extension of the dyadic solenoid, and the example of this section is an almost 1-1 
extension of the example of section [3l 

Before preceding with an analysis of As, we review some facts about the period- 
doubling substitution. This substitution forces the border |Kel] on the right, since every 
substituted letter begins with a 0. To construct the Anderson-Putnam complex [APJ, 
we need only collar on the left. The complex (call it P, for period-doubling) is shown 
in figure HI There are three collared tiles, which we call a, (3 and 7. The tile a is 1, 
preceded by a 0, the tile P is 0, preceded by a 0, and the tile 7 is 0, preceded by a 1. 
Viewed as a map on collared tiles, the substitution sends a — > 7/?, (3 — > 7a, 7 — > Pa, 
and interchanges the two vertices of the complex. 

As before, let K be the Anderson-Putnam complex for the horizontal substitution 
space X a . The approximant L for contains a piece K x K x K x [0, 1] for each edge 
a, P, 7 of P, with the identifications 

(20) (a, x, y, z, 0) ~ (J3, x, y, z, 0) ~ (J3, w, x, y, 1) ~ (7, w, x, y, 1) 
and 

(21) (-f,x,y,z,0) ~ (a,w,x,y, 1). 

The 5-tuple (a, x, y, z, t) (resp. (P,x,y,z,t) or (7, x, y, z, t)) means that the origin is 
sitting at height t in a row that corresponds to a (resp. P or 7) in the period-doubling 
substitution. Within the row containing the origin, the horizontal position of the origin 
corresponds to y G K. In the rows below and above that, the horizontal positions of the 
points one unit below and above the origin correspond to x G K and z G K. 
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The substitution E induces the following map on L, also denoted E: 



(P, (7 2 (y),a 1 (y), a 2 (z), 2t - 1); if t > 1/2, 
( 1 ,a 1 (x),a 2 (y),a 1 (y),2t); if t < 1/2. 



(22) 

E(f3,x,y,z,t) 

(23) 

S(7,a;,2/,«,t) 

(24) 



{/3, a(y) + A, a{y), a{z) + A, 2t- 1); if t > 1/2, 
( 7 ,a(a;),cT(y) + A,(T( 2 /),2t); if i < 1/2. 

(a, a 2 (i/), 02(2), 2t - 1); if t > 1/2, 

(7,<Ti(a;), a- 2 (3/), ^(j/)^*); if t < 1/2. 

(a, a(y) + A, a(y), a{z) + A, 2t - 1); if t > 1/2, 
(7> o(a;), cr(y) + A, cr(2/), 2i); if t < 1/2. 

(a, a 2 (i/), mfo), 0-2^), 2t - 1); if t > 1/2, 
(^^(a;), 0-3(2/), ai(y),2t); if t < 1/2. 

(a, 0(2/) + A, a(y), a{z) + A, 2i — 1); if t > 1/2, 
(Aa(x),a(2/) + A,a(2/),2t); if t < 1/2. 



As before, E is homotopic to a map E' that does not contain a translation by A, and 
E' = Ej o E 2 = E 2 o Ei, where E x implements a horizontally, and E 2 implements the 
period-doubling substitution vertically. Specifically, 

(25) Ei(C, x, ?/, z, t) = ((, <t(x), cr(y), <r(z),t), 
where £ = a, /3 or 7, and 

(26) E 2 (a,x,2/,2;,t) = 

(27) E 2 (/?,x,2/ ) z,t) = 

(28) E 2 (7,2,2/,^t) = 

We adopt the same strategy as in section [3l First we compute the cohomology of L 
using Mayer- Vietoris, then we take the direct limit under E*, and then we take the direct 
limit under Eg. 

To compute H*(L), let V consist of neighborhoods of the two vertices in P, and let U 
be a slightly thickened complement to V. Now U retracts to three copies of K x K x K, 
one for each edge of P, and V retracts to two copies of K x K, one for each vertex of 
P. There are 6 copies of K x K x K in the retraction of U l~l V, one at the beginning of 



(P,y,y,z,2t- 1); 


if t > 1/2, 


(j,x,y,y,2t); 


if t < 1/2. 


(a,y,y,z, 2t - 1); 


if < > 1/2, 


(l,x,y,y,2t); 


if t < 1/2. 


(a,y,y,z, 2*— 1); 


if t > 1/2, 


(P,x,y,y,2t); 


if t < 1/2. 
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each edge and one at the end of each edge. We then have: 
H°(U) = Z a ©Z /3 ©Z 7 , 

(29) H\U) = (fi x © fi y © fl z ) e 

eG{a,/3,7} 
eG{a,/3,7} 

R3 ( U ) = {fx®fy®fz)e- 
ee{a,/3,7} 

H°(V) = z 2 , 

(30) H\V) = (fl 

# 2 (\/) = (fl xy © /2 WZ ) © (fx xy © /2 yj? ) 
i/ 3 (V) = 0, 

where the two copies in H k (V) correspond to the two vertices in P. Since there is a 
copy of H k (U) at the beginning and at the end of each edge of H k (U fl V), we have that 
H k {U n V) = (H k {U)f. 

Writing down the generators of H k (U fl V) in the same order as the corresponding 
generators of H k (U) and H k (V) and keeping the copies of H k (U) corresponding to the 
outgoing edges separate from those corresponding to the incoming ones, we can again 
compute matrices for v in each dimension. The matrices will have two vertically aligned 
identity matrices in the U columns and then linearly independent vectors with two —Is 
in them in the V columns, so the signed restriction maps H k (U) © H k (V) — > H k (U fl V) 
are injective except in dimension zero. This implies that the restriction maps H k (L) — > 
H k (U) ®H k (V) are all zero (except in dimension zero), and that H k+l (L) is the cokernel 
of the signed restriction map H k (U) © H k (V) — > H k (U HV). It is simple linear algebra 
to use the generators of H k (U fl V) and the image of v to calculate that 

H°(L) = Z, 

H\L) = Z 2 , 

(31) H\L) = /i 5 , 

H\L) = (/i©/i) 7 , 

H 4 (L) = (/i©/i©/i) 3 , 

with the following generators, subject to the following constraints. 

Of course H°(L) is generated by 1. The generators of H 1 ^) are dt a , dtp and <it 7 , with 
dt a = dt y . (In other words, H 1 ^) = iJ 1 (P).) The products of generators of fl x , jX y and 
fi z with dt a , dtp and dt^ generate H 2 (L). Using the convention that generators ranging 
through ft are referred to as simply ft, we see the generators of H 2 (L) are subject to the 
four constraints 

(32) jiy U dt a — jl x U <it 7 , fi z U dt a = jx y U <it 7 , 

jiy U (dtp + dt^) — jl x U (dt a + dtp), p, z U (dtp + dt^) = fi y U (dt a + dtp). 
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The products of (p, x U jly), (jl x U jl z ) and (jly U jl z ) with dt a , dtp and di 7 generate H 3 (L), 
subject to the two constraints 

(33) fl x U ji y U di 7 = /2 y U fi z U <it a , /i x U fi y U (<it a + Gft/3) = fl y U fl z U {dtp + <it 7 ). 

Finally, H 4 (L) is generated by /i^ U /2 y U /i z U (g^ q , dig, and dt 7 ), with no constraints. 

Note the form of the constraints on H 2 and H 3 . For H 2 , we have two constraints for 
each vertex of P. The sum of the jt y U dt terms from the edges flowing into the vertex 
equals the sum of the fl x U dt terms from the edges flowing out of the vertex, and the sum 
of the jx z U dt terms from the edges flowing into the vertex equals the sum of the jl y U dt 
terms from the edges flowing out of the vertex. These may be treated as constraints 
among the fl x U dt and jl z U dt generators, while the jl y U dt generators are unconstrained. 
For if 3 we have one constraint per vertex, namely that the sum of the jlyU fl z U dt terms 
from the edges flowing into the vertex equals the sum of the fi x U jl y U dt terms from the 
edges flowing out. 

As in section [31 the direct limit of H*(L) under S* takes the same form as H*(L), only 
with fl replaced by fi. What remains is to take the direct limit under which we do 
one dimension at a time. 

The computation in dimension is trivial, and we of course have H°(X^) = Z. 

The computation in dimension 1 does not involve \x at all, and is identical to the com- 
putation of H 1 of the period-doubling substitution space. The answer is that = 
Z[1/2]©Z. 

For dimension 2, we first look at the \i y U dt terms. Since E?; maps these terms to 
themselves, 

YT 2 {jJ,yU dt a ) = fly U dtfi + fiy U rfty , 

(34) T^ifJ-yU dtp) = fiy U dt a + fiy U di 7 , 

T.Kfiy U dt 7 ) = fiy U dt a + fiy U dtp, 

(° 

these terms give the direct limit of Z 3 under the matrix I 1 

V 

eigenvalues 2,-1 and -1, and the direct limit is Z[l/2] © Z 2 . 
Next we look at the fi x U dt terms. We compute 

^(fMx U dt a ) = fly U dtp + fly U dtj, 

(35) J^KfixUdtp) = fiy U dt a + fi x U dt*y = 2fiy U dtp, 

Y? 2 (fL x U c/t 7 ) = fi x U dt a + fi x U dtp = fi y U dtp + fi y U dt^. 

Since the fi x U dt terms map to the fi y U dt terms, they do not contribute anything 
additional to the direct limit. Likewise, the ft z U dt terms map to the fi y U dt terms, and 
do not give any additional contributions. 

In dimension 3, the image of £2 consists entirely of fi x U fi y U dt and fi y U fi z U dt terms, 
since the image of £2 never involves both x and z. Moreover, the fi x U fi y U dt terms map 
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to themselves, and the jj, y U [i z U dt terms map to themselves. Specifically, we have 

U ii y U dt a ) = 0, 
E^/i^U fiyU dtp) = /j, x U ii y U cft 7 , 

(36) ^(a^x U /i^ U dt 7 ) = [j, x U fiy U (dt a + dtp), 

T^^jlyU jl Z U dt a ) = fly U fl Z U (dtf3 + dtry), 

E^/ij, U dtp) = fi y Ufi z Udt a , 
T.Kfiy U fi z U dt 7 ) = 0. 

Before taking the constraints into account, the action of £2 has eigenvalues 1, 1, —1, —1, 
and 0, making the direct limit of this action Z 4 . Because we know that Jl x U fi y U dt 7 = 
f~i y Uf~i z U dt a and p, x U fl y U (dt a + dtp) = fi y U fi z U (dt^ + dt 7 ) in H 3 (L), and because 
£2 annihilates the /i x U fi y U rft a and fi y U fi z U dt^ terms, there are only two independent 
generators left in the direct limit for H 3 (Xy,). These are [i x U fi y U <it 7 = fi y U fi z Ut a and 
^ U /ij, U dtp = jiyU ji z U dtp. The cohomology has registered the two possible fault lines, 
one between the a and 7 rows, and one between two (3 rows. 

In dimension 4, £2 is identically zero. 

To summarize, the cohomology of X s is 

H° = Z, 

H 1 = Z[l/2]©Z = iJ 1 (X pd ), 

(37) H 2 = n@ii@fji[l/2]=iJi®(H 1 (X pd )@Z), 
H 3 = (/i® /i) © (/i® fj), 

H k = for k > 3. 

Note that if 2 (X s ) is not the tensor product of fj, with the first cohomology of the 
period-doubling substitution. Rather, it is the tensor product of /i with the direct limit 
of the transpose of the period-doubling substitution matrix (as applied to collared tiles), 
which has an additional factor of Z. The two copies of /i © n in H 3 refer to the two types 
of fault lines that can occur in X s . One has a (3 row both above and below the fault line. 
The other has an a row below the fault line and a 7 row above. £ maps each of these 
situations to the other. 

5. General direct product variations with regular fault lines. 

The results of the previous section are suggestive of how the fault lines affect the 
cohomology of a tiling space. In this section we prove two theorems on the cohomology 
of tiling spaces with horizontal fault lines that arise as direct product variations. 

Suppose we have a collection of primitive 1-dimensional substitutions 07, a 2 ,---,cr N 
defined on the same alphabet A. Moreover, assume that for each a £ A, the length 
of 07(a) is the same for all k, that all of the substitutions have the same stretching 
factor, and that they all yield the same tiling space. (E.g., all of the o^s might be cyclic 
permutations of one another). Now suppose we have another primitive 1-dimensional 
substitution p. We could then consider the direct product of the two substitutions, with 
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o\ acting horizontally and p acting vertically. We then replace o~\ with etc in various 
rows, so as to introduce fault lines. Furthermore, we assume that regular fault lines occur 
at every boundary of infinite-order p-supertiles. 

We call the resulting 2-dimensional substitution S, and compute the cohomology of 
X-£. Let p = H 1 (X a ) and let v = H 1 (X p ). Let n be the number of configurations 
in the X p tiling space in which two infinite-order p-supertiles meet at the origin. By 
assumption, this is the same as the number of ways that a regular fault line can occur in 
an X s tiling. 

Theorem 1. Under these circumstances, the cohomology of is as follows: 

#°(A S ) = Z, 

H\Xv) = u, 
(38) H 2 (Xx) = p® (//©Z™- 1 ), 

H\X^) = {n®n) n , 

H k (Xx) = 0fork>3. 

Proof: As usual, let P be the AP complex of X p , and let K be the AP complex of X a . 
We take L to be one copy of K x K x K x [0, 1] for each edge in P, with the following 
identification. If the end of edge a meets the beginning of edge /3 at a vertex in P, then 
(a,w,x,y,l) ~ (J3,x,y,z,Q). 

We compute H*(L) by Mayer- Vietoris. Let V be a union of neighborhoods of the 
vertices of P, and let U contain the edges. U retracts to a number of copies of K x K x K 
(one per edge), while V retracts to a number of copies of K x K (one per vertex). UClV 
retracts to a number of copies of K x K x K (two per edge, one at the beginning and 
one at the end). 

Using the same bookkeeping as in section HI one can see that the signed restriction 
maps H k (U)®H k (V) — > H k (UnV) are all injective, except in dimension 0. The cokernel 
is understood as follows. The image of the restriction of H k (U) to H k (U fl V) merely 
identifies (up to sign) the two K x K x K terms in U n V that correspond to each edge. 
This gives obvious generators for H k (L): In dimension 1 we have dt^, where ( is an edge 
in P, in dimension 2 we have fx x U dt^, \x y U dt^, and fi z U dt^, in dimension 3 we have 
fx x U fiy U dt{ , fi x U fi z U dt^, and fi y Ufi z U dt^, and in dimension 4 we have fi x U fi y U fi z U rft^, 
all subject to the identifications imposed by the image of H k (V) in H k {U D V). In H 1 
this is that the sum of the dt^s entering a vertex equals the sum of the dt^s coming out. 
In H 2 it is that the sum of the fi y U dt terms from the edges flowing into the vertex equals 
the sum of the fi x U dt terms from the edges flowing out of the vertex, and the sum of the 
fi z U dt terms from the edges flowing into the vertex equals the sum of the fi y U dt terms 
from the edges flowing out of the vertex. For if 3 , the the sum of the fi y U fi z U dt terms 
from the edges flowing into a vertex equals the sum of the fi x U fi y U dt terms from the 
edges flowing out. For if 4 , there are no constraints. This completes the computation of 
H k (L). 

Now we note that £ is homotopic to the product of a in the horizontal direction and p 
in the vertical direction. Taking the direct limit under a* merely converts fx to p. What 
is left is taking the direct limit under p* . 
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In dimension 1, this gives v. 

In dimension 2, we note that the constraints express certain combinations of the p x U 
dt^s or the p z U dt^s in terms of the p y U dt^s. They do not constrain the p y U dt^s terms, 
which are mapped to themselves by p*. Furthermore, among these terms the action of 
p* is just the transpose of the substitution matrix itself (as applied to edges of P, i.e. to 
collared tiles). The direct limit of H 2 (L) under p* therefore contains the direct limit of 
this matrix. 

In principle, the direct limit of H 2 (L) should also contain contributions from the p x Udt^ 
and p z U dt^ terms. However, we claim that these contributions are zero as a consequence 
of our using a complex that forces the border, as explained below. 

Note that the pullback of dx U dt^ term is a dx U dt$ term for each tile £ for which 
is a word beginning with (, plus a dyUdt^ term for each tile £ for which contains ( in 
the middle or end of the word. If the substitution forces the border in m steps, then p m 
of each edge emerging from a vertex in P is a word beginning with the same letter (call it 
a;). (p m )*(p x U dtj) then equals the sum of the p x U dt terms from all the edges emerging 
from this vertex (plus additional p y U dt terms). However, the sum of the p x U dt terms 
from the edges emerging from the vertex equals the sum of the p y U dt terms from the 
edges entering the vertex. Likewise, if 77 7^ u>, then (p m )*(p x U dt v ) contains none of the 
p x U dt terms from edges emerging from the vertex. Either way, the pullback by p m of a 
p x U dt term can be expressed as a sum of p y U dt terms. The same argument (applied 
to ends of words) shows that the pullback by p m of a p z U dt term can be expressed as 
a sum of p y U dt terms. In particular, the p x U dt and p z U dt terms that are linearly 
independent of the p y U dt terms do not appear in the eventual range of p*, and hence 
do not contribute to the direct limit of H 2 (L). 

In dimension 3, p* sends p x U p z U dt terms to p x U p y U dt and p y Up z U dt terms. 
Therefore, we need only consider the direct limit of p x U p y U dt^ and p y U p z U dt^ terms. 
We associate p x Up y Udt^ with the vertex in P that ( leads into, and associate p y Up z Udt^ 
with the vertex it leads out of. As in dimension 2, forcing the border implies that (p m )* 
takes each p x U p y U dt term to either all or none of the p x U p y U dt terms associated 
with a vertex. When dealing with p x U p y U dt, we may therefore restrict our attention 
to sums of all the edges emerging from a vertex, and when dealing with p y Up z U dt we 
may restrict our attention to sums of all the edges leading into a vertex. However, the 
sum of all the p y U p z U dts from the edges leading into a vertex equals the sum of all 
the p x U p y U dts leading out of the vertex, so we have exactly one independent term per 
vertex. 

The substitution p maps vertices of P to vertices, and has a natural pullback action 
on the 3-forms associated to vertices. Eventually, p merely permutes the vertices that 
describe boundaries of two infinite-order supertiles. The 3-forms associated with these 
vertices give a basis for H 3 (Xj^). 

In dimension 4, the pullback map is zero, so if 4 (Xs) = 0. 

Finally, we consider the coefficient of p in H 2 (X-s) and the coefficient of p <S> p in 
H ?> {Xy). The first is the direct limit of the transpose of the substitution matrix as applied 
to edges of P, and the second is the direct limit of the transpose of the substitution matrix 
as applied to vertices of P. In other words, they are the direct limit of 1-cochains and 
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Figure 5. The Anderson-Putnam complex for the vertical tiling space 

0-cochains on P under substitution. These are closely related to direct limits of the 
cohomology of P (i.e., to the cohomology of X p ). Since the direct limit of the 0-cochains 
(namely Z n ) is Z n_1 more than H°(X P ), the direct limit of the 1-cochains must be Z n_1 
more than H 1 (X p ). | 

In stating Theorem 1, we assumed that every boundary between infinite- order super- 
tiles in X p led to a fault line. As the following example shows, this is not always the 
case. 

Let Ai, A 2 and A% be tiles of width A and height 1/3, and let Bi, B 2 , and B 3 be tiles 
of width 3 and height 1/3. Our substitution is 

(39) E(A) = (* *) ^ = {t *) E( - 4 » ) = {t ft) 

Each row in a tiling contains either Ajs and i?is (call this an a row), v4 2 s and -B2S or 
/I3S and B 3 s. The complex P consists of three edges (a, j3, 7) running in a circle, with 
a followed by j3 followed by 7 followed by a. The vertical substition p is 

(40) p(a) = a/3, p(J3) = 7a, p{l) = h- 

Fault lines do not develop at the boundary of a and (3 supertiles, or at the boundary of 
(3 and 7 supertiles, just at the boundary of 7 and a. 

This example may look mysterious, but it is just a rewriting of the basic example of 
sectional Ai, A 2 , and A 3 are just the bottom, middle and top thirds of the A tile, while 
Bi, B 2 and B 3 are the bottom, middle and top thirds of the B tile. The cohomology is 
correctly described by Theorem 1, only with n being the number of possible fault lines 
(namely 1), not the number of ways that two infinite-order supertiles can meet (namely 
2). This observation generalizes to 

Theorem 2. Suppose we have a 2- dimensional substitution E ; generated by a vertical 
substitution p and horizontal substitutions o~\, . . . , o~n, as in Theorem 1. Suppose that 
the boundaries between infinite- order p- supertiles are either fault lines or are rigid, with 
the patterns on both sides of the boundary being mutually locally derivable. As before, 
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let p = H 1 (X cr ) and let v = H 1 {X P ). Let n be the number of ways that a fault line can 
develop in an X-% tiling. Then the cohomology of X^ is as follows: 

ff°(A s ) = Z, 

H\Xv) = 

(41) H 2 (Xx) = p® (//©Z"- 1 ), 

ff 3 (X s ) = {p®p) n , 

H k (X s ) = 0fork>3. 

Proof: Let P be the Anderson-Putnam complex of p. By assumption, each vertex of P 
either generates a fault line (call this an essential vertex) or has the patterns on both 
sides of the vertex precisely aligned. We rewrite p using the essential vertices as stopping 
and starting rules as in |BDlj . and rewrite £ in terms of these new vertical tiles. By 
construction, each of the n vertices of the new vertical substitution generates a fault line, 
so Theorem 1 applies directly. ■ 

6. Open problems 

(1) We understand fault lines for substitutions on two letters, but what about more 
complicated substitutions? Is it possible to have lines without finite local complex- 
ity that do not allow arbitrary shears? What can we say about the cohomology 
of tiling spaces that allow such "irregular" fault lines? 

(2) In considering Theorem [2j we assumed that all boundaries between infinite-order 
p-supertiles either generated fault lines or had the two sides remain in lockstep. 
When the horizontal stretching factor is not Pisot, are these the only possibilities? 

(3) What happens if the horizontal substitutions of Theorem 1 are different enough 
that the rows in a supertile are not all the same up to translation? (The hori- 
zontal stretching factors would all have to be the same, which would constrain 
the possible abelianizations of the different cxj, but the actual substitutions could 
differ.) There are natural conjectures for what if 1 and if 3 of such a tiling space 
should look like. if 1 should come entirely from the vertical substitution and if 3 
should contain a copy of p± ® pi for each possible fault line, where p\ (resp. p^) 
is if 1 of the tiling space that describes the row immediately above (below) the 
fault line. However, it is not at all clear what if 2 should be. 

(4) Some tilings (for instance the one in Figure 2) allow both vertical and horizontal 
fault lines. Since a single tiling can exhibit a lack of coordination across at most 
one fault line, it is easy to guess what if 3 of such a tiling space should look like, 
with a contribution (as in the previous problem) from each possible fault line, 
vertical or horizontal. But does if 1 vanish altogether? What about if 2 ? 

(5) Up to this point we have only been considering rectangular tiles. What if the 
tiles are not rectangular, as in the generalized pinwheel? Again, it is possible to 
predict the cohomology in the highest dimension (4 for the generalized pinwheel, 
since the rotations provide an additional degree of freedom), with a contribution 
from each possible species of fault line. However, we do not venture to guess the 
lower dimensional cohomology. 
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